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LONGITUDINAL RESISTIVITY AND INTERMEDIATE MACH NUMBER 
COLLISIONLESS TRANSVERSE MAGNETO SO NIC SHOCKS 

I. Introduction 

It is now reasonably well established that transverse collisionless 

shocks can exist in at least two Mach number regimes. For m (the 

magnetosonic Mach number) less than about 2.5, transverse resistivity is 

the dissipation mechanism which forms the shock.Since the plasma is 

nearly collisionless, this resistivity is anomalous in that it is assumed 

to arise from an electron-ion streaming instability in the shock front. 

Measurements of the turbulent spectrum in the shock front show that 

fluctuations with wavenumber k ~ kp, the Debye wavenumber, are indeed 

excited. * It is generally assumed that the resistive dissipation heats 

only the electrons. However a small amount of ion heating cannot be ruled 

out. It is certainly not easy to see how an instability can heat only one 

species. An ion acoustic instability, at least, does give rise to a small 

amount of ion heating.^ In fact one theory of the resistive shock 
a 

structure has this small ion heating playing an important role in 
determining shock structure. 

For larger m, it is now well established both from theory^*^ and 
12—19 

particle simulations, 4 that ion reflection from the shock front can act 
as a dissipation mechanism in one of two ways. On the one hand, these 
reflected ions may stream ahead of the shock if they are unmagnetized or 
sufficiently weakly magnetized.^*^®*12-15 if this be the case, the 
standard Rankine-Hugoniot conservation relations do not apply because of 
the additional species upstream. On the other hand, the reflected ions may 
turn around in the magnetic field and drift downstream, their large Larmor 
orbit velocity effectively becoming part of the downstream 

pressure.**»16-19 Since there is no additional upstream species in this 
Manuscript approved October 27, 1983. 
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case, Rankine-Hugoniot conservation relations do apply. In fact the 
simulations of Refs. 18-19 were successful because the Rankine-Hugoniot 
relations were exploited in initializing the system. Of course for the 
latter dissipation mechanism to be significant, the time of the experiment 
or simulation must be longer than about ir/u) c j, where is the ion 
cyclotron frequency. 

The purpose of this paper is to show that between these two ranges of 
Mach number, there exists an intermediate range. In this intermediate 
range, the shock structure and dissipation are dominated by the breakdown 
of the quasi-neutral approximation. Above a critical Mach number, the 
quasi-neutral fluid equations predict infinite slope, Implying, of course, 
the breakdown of quasi-neutrality. However, the full set of fluid 
equations, including charge separation, are nonsingular. Thus, above this 
Mach number, the flow sets up an electrostatic ion acoustic wave train. 
The Landau damping of this wave is the dissipation mechanism which forms 
the shock. This Landau damping may heat either the electrons or ions. In 

purely electrostatic shocks with very cold ions, the electron Landau 

20 

damping seems to dominate. Since the electron Landau damping rate is 
small, the shock structure is oscilliatory. On the other hand, in 
magnetosonic shocks the ions are warmer, so ion Landau damping dominates 
and the shock structure is monotonic. It is interesting that this 
dissipation mechanism does not rely on exciting some plasma instability in 
the shock front, and in that sense it is not anomalous. Indeed, the 
electrostatic subshock has a Debye length scale size, so it is almost 
impossible that an instability with still smaller wavelength could be 
excited. 
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While this dissipation mechanism is not anomalous it is very 
complicated. The simplest model of it shows tnat it is nonlocal, and in 
that sense non fluid like. However, there is a new simple fluid 
dissipation model. When the electrostatic oscillation is set up, the 
longitudinal electron and ion velocities are no longer equal. A 
longitudinal friction force between the two is an additional possible 
dissipation mechanism which apparently has not been examined in the 
previous literature on high Mach shock waves. An important purpose of this 
paper is to investigate the effect of such longitudinal resistivity on 
shock structures. Our conclusion is that it is a very good fluid model for 
collisionless magnetosonic shocks with Mach number in the intermediate 
regime. It has a precise onset, namely when the shock reaches the critical 
Mach number. It shows that shocks above this Mach number have an 
electrostatic subshock embedded within in and its structure is governed by 
longitudinal resistivity. When the Mach number increases to a second 
critical value, the full fluid equations become singular, and this defines 
the upper limit of the intermediate range. The interpretation of this 
second singularity is simply made in terms of ion reflection. Above this 
point, the shock, if it exists, is governed by ion reflection as described 
in, for Instance, Refs. 9-19. At no point is it either necessary or 
physically correct to postulate a viscosity. 

While the theory we deal with here concerns only steady state shocks, 
the general conclusions should also apply to non-steady flow of 
Interpenetrating plasmas. Specifically, electrostatic interactions should 

be Important in a recent series of laboratory experiments on the flow of a 

21 

laser produced plasma through a background magnetized plasma. 
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In Section II we review experimental data on high Mach shock 

structures and show that it indicates the existence of this intermediate 

regime* To do so, we note that in this Mach number range many of the 

experiments are over in a time of order of or less than tt/oj . , so 

cl 

1 3 22 23 

convection of reflected ions downstream * »» cannot be the dissipation 
mechanism of importance* Also, downstream ion temperatures are measured in 

several experiments and they agree with Rankine-Hugoniot. This means that 

3 22 

pure reflection of unmagnetized ions, ' cannot be the dissipation 

mechanism either. 

In Section III we discuss shock structure for unmagnetized ions where 
the quasi-neutrality assumption is valid and also where it is not. We find 
that for sufficiently weakly magnetized plasmas, characteristic of the 
experiments in Refs. 1-6, 22, 23, the electrostatic and magnetic structures 
essentially decouple, and that the electrostatic subshock must ride on the 
crest of the magnetic shock structure. This also agrees with experimental 
results. 

In Section IV, and also in Appendix B, we discuss the actual 
dissipation mechanism (Landau damping) and fluid models for it. As we 
show, even the simplest approximation to it is quite complex, and does not 
easily fit into a simple fluid model. However the actual physical model 
provides Important insight on a best guess for longitudinal resistivity and 
on the partition of resistive heating between electrons and ions. Another 
interesting phenomena is that Landau damping can only be a valid long time 
description of the dissipation mechanism if the shock structure fluctuates 
in time. Thus the shock cannot be time independent but, of course, its 
average profile can be. 
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Section V presents the results of shock structure calculations in the 



intermediate Mach number regime. Also, where possible, comparisons are 
made with experimental results. Finally, Seccion VI draws conclusions. 
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II. Review of the Results of Laboratory Experiments 


on Transverse Collisionle3s Shocks 


In this Section we review briefly previous experiments of transverse 

1—22 2 ^ 

shocks in laboratory plasmas. °>*- 4 ** J Our conclusion, based on 

examination of this data, is that for m < 2.5, the shocks are resistive, as 

described by authors of these references. However, there is a 

range 2.5 < m < 5 where steady shocks form which are above the critical 

Mach number for resistive shocks. However, the process of ion reflection 

followed by downstream convection by the E x B drift cannot be operative 

either, if only because the duration of the experiment is less than the ion 

Larmor time. Furthermore, the dissipation mechanism cannot be simple ion 

reflection of magnetized ions either. For ion reflection in an 

unmagnetized plasma, the second ion species upstream drastically alters the 

Rankine-Hugoniot conditions downstream.^ However, in two of the 
o 22 

experiments cited, downstream ion temperatures were measured and found 

to agree well with Rankine-Hugoniot conditions. Thus our conclusion is 
that there must be some sort of new dissipation mechanism associated with 
this intermediate range of Mach numbers. 

In Ref. 1, Paul et al . launch shocks with Mach numbers 2.5, 3.7 and 
6.3. For the first, resistivity forms the shock and gives the electron 
heating required to satisfy the Rankine-Hugoniot relations. In the case 
of bi • 3.7, the downstream electron temperature is less than that predicted 
by Rankine-Hugoniot, implying that the ions are heated. In addition to the 
normal magnetic structure, there is a very low amplitude magnetic pedestal 
in front of the shock. This pedestal is generally interpreted to mean a 
small number of ions are reflected from the shock front and convect the 








magnetic field upstream. However in Paul's experiment, the convection of 
these ions downstream by the E_ x B_ drift is almost certainly not the 
mechanism which forms the shock, at m * 3.7. The shock velocity is about 
2.5 x 10 7 cm/sec and the shock propagates ahead of the piston a distance of 
10-15 cm, implying a time of about 400 nsec. However, ir/co^ is also about 
400 n sec. Thus any reflected ions cannot make it back behind the shock 
until the shock hits the cylindrical axis. However, shocks are observed 
for much shorter times. For the measurement at M * 6.3, the upstream 
pedestal broadens and increases in amplitude so that it encompasses the 
entire shock. At this Mach number, reflected ions are almost certainly 
playing a role. 

o 

In another series of experiments, Keilhacker et al. measured the 
properties of magnetosonic shock waves with Mach numbers between 1.5 and 
5. Here ion temperatures behind the shock are measured by collective 
scattering of laser light at wavelengths less than the Debye length. The 
measured ion temperatures agree with Rankine-Hugoniot conditions. In this 
experiment, the total transit time for the shock is about ly sec. Thus the 
dissipation mechanism cannot be the convection of reflected ions downstream 
because there is not enough time; neither can it be the pure reflection of 
unmagnetized ions, because the downstream ion heating corresponds to that 

predicted by Rankine-Hugoniot. 
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Segre and Martone have also investigated transverse shocks in the 
Mach number range 3.1 < m < 4.25. In their experiment, there is no 
magnetic field pedestal leading the shock, so reflected ions almost 
certainly play no role. Even if there were reflected ions, the ion period 
- 450 n sec, the propagation time is only about 250 n sec, so no 
reflected ions could make it back downstream. In Ref. 22 the downstream 











electron density and temperature are measured by laser scattering. However 
magnetic measurements are used to measure n(T g + T^) downstream, so that 
downstream ion temperatures can be directly inferred. The measured ion 
temperatures do in fact agree with that predicted by Rankine-Hugoniot. 

23 

Another series of shock experiments are those of Eselevich et al. 

The shocks are produced in a tube of 8 cm radius in H, He or A. Since the 

13 15 -3 

gas can be varied, and the density varies between 10 < n < 10 cm and 

100 G < B < 2000 G, there was a very large parameter space studied. 
Eselevich et al. found that for Mach numbers between 3 and 5, steady she s 
form. In this range of Mach numbers, a pedestal in magnetic field might 
might not lead the shock. For instance, Fig. 2a of Eselevich et al. si 
a Mach 4 shock in helium with no pedestal, whereas Fig. 4 shows a Mach 4 
shock in hydrogen with a small (i.e., less distinct than in Paul's 1 
experiment) magnetic foot leading the shock. Thus, reflected ions are not 
a universal, and probably not a very important, aspect of steady transverse 
shocks in this intermediate range of Mach numbers. However, Eselevich et 
al., like Paul et al., did observe that as the Mach number increases 
(beyond about 5.5 generally in Ref. 23), the magnetic foot does expand out 
and dominate the entire structure. 

The moat distinct feature of Eselevich's observations are that in the 
intermediate range of Mach numbers, an isoraagnetic potential jump forms on 
the crest of the shock wave. The thickness of the potential jump is 
typically a Debye length. Thus, a key feature in this intermediate range 
is the breakdown of the assumption of quasi-neutrality. The potential jump 
across this isomagnetic transition region initially increases with Mach 
number, but then above some higher Mach number begins to decrease until 
above a second critical Mach number, typically about 5 or 6, the entire 











isomagnetic jump disappears and the shock structure is taken over by the 
magnetic pedestal. Hence, we feel that a key element in a fluid 
description of shocks in this intermediate range is the breakdown of quasi¬ 
neutrality. Once quasi-neutrality is violated, an additional fluid 
dissipation mechanism, longitudinal resistivity becomes possible. In the 
next two sections we examine the consequence of longitudinal resistivity. 
Although the actual dissipation mechanism is considerably more complicated 
and non fluid like (see Sec. IV), a longitudinal resistivity model can 
explain many features of the experimental results. Specifically it hinges 
on the breakdown of quasi-neutrality; it provides for the proper electron 
and/or ion heating; it predicts an isomagnetic potential jump on the crest 
of the shock wave; and it gives the shock structure only in an intermediate 
range of Mach numbers between the critical Mach number for transverse 
resistivity and the onset of strong ion reflection. Thus, longitudinal 
resistivity is probably as good an approximation to the dissipation 
mechanism that one can obtain within the limits of a conventional fluid 


description 


















III. Shocks with Transverse and Longitudinal Resistivity 


In this Section, we discuss shocks formed by longitudinal resistivity. 
Before doing so, however, in Section III.A we set the stage by reviewing 
the structure of shocks with transverse resistivity. Section III.B 
discusses the effect of longitudinal resistivity, and Section III.C sets 
out the procedure for calculating the shock structure in weakly magnetized 
plasmas. 
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A. Review of Shocks Formed by Transverse Resistivity 

To start we formulate the equation for quasi-neutral transverse 
resistive shocks. We neglect both transverse and longitudinal electron 
inertia, and the motion of the ions in the transverse direction. Assuming 
that the shock propagates in the x direction, the magnetic field is in the 
z direction, and upstream parameters have a zero subscript, then the 
electron and ion density equation gives the result 


nV - n V , 
o o’ 


( 1 ) 


where a velocity or electric field without a direction subscript means a 
velocity or electric field in the x direction. Then the electron momentum 
equation in the x and y direction give respectively equations for E and 


V ey . They are 


- i— n T - B 
ne dx e c 


( 2 ) 


and 


ey 


eV n 

-2. (b-— B], 

m\jr v n J 9 


(3) 


where m is the electron mass, v the effective collision frequency, and the 
resistive force density on the electrons is - nmv(V ev - ] = - nmv V 

since we assume V 


ey 


ey iy' " "' ey 

0 (unmagnetized ions). Also we have used the fact 


that Ey - E yo - — B o . 


To write out the electron temperature equation we first assume that 
the electrons have y « 5/3. Then we must also specify what fraction of the 

















resistive heating goes into the electrons and what fraction goes into the 
ions. In terms of the ion y, the total energy equation for electrons plus 
ions is 


4 - fi nMV 3 + 1 nVT + -X- nVT. ) - - neV E , (4) 

dx ^2 2 e y-l i 1 ey y’ v ' 

whereas the total energy equations for the electrons and ions separately 
are 


nVT 

dx 2 e 


- neV »E + A 


(5a) 


and 


(1 nMV 3 + ^ nVT t ) « neVE - A, 


(5b) 


where M is the ion mass. The term A denotes an energy exchange between 
electrons and ions, and its value depends on what fraction of the resistive 
heating goes to electrons and what fraction goes to ions. To determine 
this, one must write out the temperature equations for electrons and 
ions. To do so, take the dot product of the electron momentum equation 
with Vg and subtract the result from the electron energy equation to get 


j nVT 
2 dx e 



A + nmv 


V 2 . 
ey 


(6) 


Then multiply the x component of the ion momentum equation 

















by V, subtract the results from the ion energy equation to obtain 


4 ~ nVT - V ~ nT 
2 dx i dx i 


- A. 


(8) 


If a is the ratio of ion resistive heating to electron resistive heating, 
then 


or 


o 


- A 


A + nmvV 2 
ey 


(9a) 


\ 


A 


a 

1 + a 


nmv 


V 2 . 


ey 


(9b) 


The quantity a depends on the physical process causing the resistive 
heating. Classical resistivity gives a - m/M. Anomalous resistivity gives 
an a which depends upon the particular mechanism. Generally a « 1; for 

Instance if the resistivity is caused by an ion acoustic instability, 

7 8 

o - 0.1-0.2. Finally, these equations are closed by Ampere's law 


dB m 4irne v 
dx c ey 


( 10 ) 


In order to see where the critical Mach number is, a standard 
procedure is to linearize Eqs. (1-10) about the upstream or downstream 
state. Assuming that ■ k, and an unsubscripted variable denotes either 
the upstream or downstream state, there are only two roots for k, 
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k - 0 


(11a) 



k > 0 so Che linearized solution diverges away from Che state. Equation 
(12) corresponds to the upstream state of the shock wave for which the flow 
speed must exceed the magnetosonlc velocity. The second regime corresponds 
to 




2 ^i 5 ^e 

< V +_- + — — 

A M 3 M 


(13) 


where k < 0. Since the flow speed is less than the magnetosonlc speed, it 
corresponds to the downstream state, and since k < 0, perturbations 
converge into the downstream state. This then is a standard (transverse) 
resistive shock having a subcritical Mach number. The final regime 
corresponds to 
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having k > 0. Since perturbations now diverge away from the downstream 
state, there can be no shock in the regime defined by Eq. (14) if 
transverse resistivity is the only dissipation mechanism. This then 
defines the critical point above which no shock can form. 

To continue, let us write out the equations for a quasi-neutral shock 
as a system of nonlinear first order differential equations, as is 
appropriate for standard numerical integration. 

The equations for B and V g y are already expressed in the proper format 
in Eqs. (10) and (3). Equations (6a) and (6b) can easily be reduced to 


dT 

e 

dx 



(16) 


and 


dT. i * T, , 

_i _ 1 r A _i dll 

dx y- 1 n V n dx 
' o o 


(15) 


but now dn/dx appears on the right hand side. To solve for dn/dx, use the 

ion momentum equation, insert the electric field from Eq. (2), use Eqs. 
dT E dT i dn 

(15) and (16) for — and -s— , and solve for -=— . 


2 2 
_ n mvV 
2 _ey 


ne —— B + -=•-- 

c 3 n V 


The result is 






















This can be used on the right hand side of Eqs. (15) and (16). Notice 

d dT i dT 2 

that —and E all become infinite at V - (yT^ + -j T g )/M. Thus, 

2 

quasi-neutral shocks formed by transverse resistivity only exist if V < 
(yT^ + y T^j/M everywhere. If the shock is monotonic and y * y (or more 
generally, if y g * Y^)» this critical Mach number can be determined by the 
Rankine-Hugoniot conditions only. However, if the shock is not monotonic, 
and/or y g * Y^» the critical Mach number depends on the details of the 
oscillatory structure and how the resistive heating is partitioned between 
electrons and ions. 

A key question then, for the description of shocks in the intermediate 


regime, is the physical significance of the singularity. One might think 

that the singularity indicates the onset of ion reflection, but this is not 

the case. To see this, imagine that T^ << T g . Then ion reflection occurs 
1/2 

where V ~ (T,/M) - 0. However, the singularity occurs at a much higher 

5T e 1/2 

velocity, V • ) . The significance of the singularity lies not in 

ion reflection, but rather in the breakdown of the quasi-neutral 

dE 

assumption. This is clear from the fact that both E and ^ , as expressed 
in Eq. (3), become infinite at the singularity, implying a breakdown of 
quas1-neutrality. Thus, in the intermediate regime, charge separation 
begins to play an Important role in the shock structure. 


















B. Shocks with Longitudinal Resistivity 

In Chis subsection and the next, we discuss how shocks are formed by 
longitudinal resistivity. If quasi-neutrality is no longer assumed, the 
steady state ion fluid equations are 


and 


n.V. - n V 
i i oo 


h («X + °l T i) ' V E • - V e>' 


h (7 “i m i + Vi T i) ■ Vi eE + V '* 2 


ey 


(18) 


(19) 


+ TTt 1 W W • 

l 


( 20 ) 


where and are defined analogously to a and v. Notice that we allow 
different values for longitudinal and transverse dissipation and energy 
partition. This is reasonable because, as we will see, these two types of 
dissipation arise from completely different physical processes. 

The fluid equations for the electrons are 


n V - n V , 
e e 00 


A v B 

j — nT - -neE-ne e ^ - + n.mv.fv. - V ), 
dxee e e c iJl'i e' 


( 21 ) 


( 22 ) 






















at k » 0, what remains in Eq. (27) is a cubic equation for k. As we will 

2 2 

see, one of the roots is of order u» V/vc , while the other two roots are 

pe 

2 1/2 2 2 
of order k D = (4im e /T e ) . Assuming that k Q » co^V/vc , the equation 

separates. The low k root is still given by Eq. (lib). We denote this the 


Magnetic Mode. To get a simple form for the high k roots, we further 
2 2 

assume that moj / Mk^v V « 1, and u> /v << 1. Since the length of a 
c u o c Jt 

resistive shock* - ^ is typically 10c/o)p e » an approximate order of magnitude 

f kk. ,. .. in.. IT / „ Tk« 
















v V _ » 

4- (4 t . + ~ mv j 


t - ^ T i - 


10 

3“ 


l _ 

(yT^V 2 ) 


K^[| T e + 1* t ' MV oD 2 + 


l 2 -^(y T e + ^ - W 2 ) [y\ - MV 2 )}' 72 . (28) 

where u g is the electron thermal velocity u g * (T /m)^ 2 . 

Notice that the wave number given by Eq. (29) is of order k D . Also, 
the magnetic field completely decouples from this large k mode, so it is 
entirely electrostatic in nature. Depending on the values of velocity and 
temperature, the linearized modes structure themselves as shown in Table I. 

Now consider how to integrate Eqs. (18)—(26) to solve for the shock 
structures. Let us say we start upstream and integrate toward the 
downstream state. Referring to Table I the upstream plasma is in Class A, 
since the flow speed must exceed the magnetosonic speed. Hence there are 
two linearized perturbations which increase away from the upstream state. 
The relative amplitude between the two is arbitrary at this point. Now let 
us assume the downstream plasma is in Class B or C. In each case there are 
two linearized eigenfunctions which converge into the downstream state, and 
one which diverges. Thus the relative amplitude of the two linearized 
perturbations upstream must be chosen so as to eliminate the divergent 
perturbation downstream. This is then a nonlinear eigenvalue problem, and 
this concept has been applied to both conventional shocks 2 ^ and reflected 
ion shocks. 
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Table I 


Velocity Regimes and linearized modes in each regime. 


Class 


Velocity Range 


Modes 




MV? + 4- T + vT, < MV 2 
A 3 e i 


One electrostatic mode 

and 

the 

magnetic 

mode 

diverge 

from 

the 

state as 

X 

increases 

in 

the 

direction 

of 

flow. 

The 

other 


electrostatic mode converges. 


B j T g + yT i < MV 2 < 

'"a + T T e + * T i 

- ■ 

2 5 

C yT^ < MV < y T g + 3T^ Magnetic mode diverges, both 

electrostatic modes converge. m 


Magnetic mode converges, one 
electrostatic mode converges, one 
diverges. 


MV < yT. 






Magnetic mode and one electric mode 
diverge, the other electrostatic 
mode converges. 
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If the downstream state is Class B, the scheme is trivial. The 
divergent downstream perturbation is electrostatic so the solution is just 
to eliminate electrostatics and calculate shock structure based on quasi¬ 
neutrality. If the downstream state is Class C, the problem is more 
complex. Here the divergent state is magnetic and both convergent states 
are electrostatic, so there is necessarily a coupling between the magnetic 
and electric structure somewhere in the shock. In the next subsection we 
will show how to calculate this coupled structure in a simple way by 
exploiting the fact that the length scale for the electric structure is 
very small compared to the length scale for magnetic structure. 

Finally, if the downstream state is Class D, there are two divergent 
states downstream, but only one free parameter upstream, so that 
longitudinal resistivity does not allow a shock to form in this case. In 
the next subsection we will discuss the meaning of the singular behavior in 















Calculating the Double Structured Shock 


As we have seen in the previous two subsections, when the plasma is 
making the transition from Class B to Class C, the electrostatic 
approximation breaks down and large electric fields and charge separations 
are set up by the flow. If the plasma were in Class C, and there were no 
dissipation, this charge separation results in an undamped ion acoustic 
wave. If it were in Class B with no dissipation the charge separation 
manifests itself as an ion acoustic soliton. With dissipation in the form 
of longitudinal resistivity, an ion acoustic shock can form so that the 
plasma can make a transition from Class B to Class C. Thus near this 
point, a narrow ion acoustic shock is embedded in the broader magnetic 
structure. This then corresponds to the isomagnetic potential jump 
discussed in Ref. 23. 

The choice of the electrostatic wave amplitude upstream is then 
equivalent to a choice of the position of the ion acoustic subshock within 
the broader magnetic structure. As we have seen, for a Class C plasma 

there is no magnetic perturbation which converges into the downstream side 

2 2 

of the ion acoustic subshock. This means for the case that k n >> oj V/vc , 

u pe 

where the magnetic structure and ion acoustic wave decouple, the ion 
acoustic wave must be at the position where the magnetic field has its 
downstream value. 

Hence the following scheme for computing the shock structure suggests 
itself. First integrate the equations for a quasi-neutral resistive shock, 
as written in Section III.A until one reaches the position where the 
magnetic field has the downstream value predicted by the Rankine-Hugoniot 
conditions. We find that at this point, the densities and temperatures are 
always below the downstream values, and also MV > 5/3 T g + yT^. At this 



















point simply set the magnetic field equal to its constant value, set V g y * 
0 and solve for the ion acoustic shock in the unmagnetized plasma. The 
equations of the ion-acoustic shock are 
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The ion density and momentum equation are as in Eqs. (18) and (19). For 
the electrons 
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and Poisson's Equation is still as given in Eq. (26). At the onset of the 
ion acoustic shock, the electric field can be initialized by taking its 
value from the quasi-neutral solution at the switch point. Of course, the 
exact value that E is initialized at is not very important because E 
reaches much higher values in the center of the transition region. 

Then Eqs. (18, 19, 26, 29-32) are integrated from the start of the ion 
acoustic shock on through to the downstream state. It is easy to show that 
this procedure preserves the overall Rankine-Hugoniot conditions; this is 
shown in Appendix A. 

To proceed, let us examine the singular nature of these equations by 

dT i 

writing out equations for V^. Eliminating from the energy equation 


(Eq. (29)) and inserting it in the momentum equation, Eq. (19), we find 




















( 33 ) 




Now the only singularity is at * yT^/M. It is especially easy to draw 

insight from this for the case y ■ 3, one dimensional ions. As shown in 

Ref. 10, the fluid equations are those derived from a waterbag ion 

1/2 

distribution for which the thermal width is (3T /M) . Then the 

2 

singularity at V * 31^/M simply means that an ion at the outer edge of the 
waterbag distribution function has been reflected (its velocity has gone to 

zero in the shock frame). The obvious generalization then is to regard the 
2 

singularity at V * yT^/M as the generalization of the condition for ion 
reflection to two or three dimensional ion distributions. Hence this 
singularity does have a simple interpretation in terms of reflected ions. 
In Ref. 10, it was shown that for ion acoustic shocks, one could add a 
reflected fluid upstream. The density of this reflected fluid then becomes 
an additional upstream parameter, which can be varied so as to eliminate 
another divergent downstream solution. This might also be possible to do 
for the magnetosonic shocks which we consider, but it will not be examined 
further here. 

To summarize, we have shown that longitudinal resistivity can give 
rise to double structured shocks in an intermediate range of Mach number. 
These are where the downstream state is Class C, as defined by Table 1. At 
the upper limit of this range, the shock structure once again becomes 


singular, but this singularity can now be interpreted in terms of reflected 
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In the previous section, we discussed how longitudinal resistivity 

forms the transverse shock in the intermediate Mach number regime (Class 

C). The idea is that as the Mach number increases, large charge separation 

is generated, so an ion acoustic wave is set up. As this wave damps by 

longitudinal resistivity, the shock forms, as discussed in Section II. 

However, longitudinal resistivity is only a fluid model for the 

dissipation, the actual dissipation mechanism is electron or ion Landau 

damping, which is inherently non-fluid like. Nevertheless, it is an 

extremely powerful dissipative mechanism. If the ions are cold the damping 

distance due to electron Landau damping is extremely small, of order 

(M/m)*^ k This implies a v of order oi • If the ions are warm, the 
D i pe 

Landau damping rate is larger, implying still larger values of 

As the longitudinal resistivity appears first in the momentum 
equation, we proceed by calculating the term in the electron momentum 
equation which arises from electron Landau damping. This means getting the 
correction to the pressure tensor. To do so we assume the plasma is 
infinite and homogeneous work in the (u>, k) domain, and then write a 
linearized fluid equation. The only term which is not completely 

straightforward is the pressure tensor. As usual, we denote perturbed 
quantities with a superscript tilda. Then 
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where now a v denotes a particle velocity in a Vlasov model and, as usual, 
a capital V is the fluid velocity. The total perturbed electron density is 
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where we have assumed Maxwellian electrons, f g ■ (/2ir u^) exp-(v /2u e ) and 
we work now in the plasma rest frame. Taking the resonant part of the 
second term in the brackets of Eq. (36), that is 
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we find 
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To continue, we would like to find the perturbed part of the total 
electron momentum flux 


P g - m / dv v* f e - - / dv 
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(39) 


The resonant part Is small by a factor of fr-^—) and we neglect it, so that 

























To obtain an equation of state, we need the perturbed pressure in terms of 
the perturbed density. From Eqs. (38) and (40) we find 



(41) 


The force on the electrons is -ikt^. The first term in the square brackets 
simply gives - T g — n g , the pressure gradient force on an isothermal 
plasma. The second term, which is out of phase by ir/2 is the additional 
dissipative contribution to the electron force. Using the fact that the 
linearized density equation gives the result n^ - u/k n e , we find the 
additional force density on the electrons is 
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Modeling this additional retarding force as a longitudinal resistivity in 
the subshock, we find 
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where we have assumed that within the subshock k . kp. 

Although the dissipation from Landau damping is quite simple in the k 
domain, it is very complicated and nonlocal in the x domain. In the k 
domain, the dissipative part of the perturbed pressure is the product of 
-i nm u g -|-£-p 1^ the x domain, it is the convolution of the 

Fourier transforms of the two factors. Thus the dissipative part of the 
perturbed pressure in the x domain is 
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where P denotes the principal value. Thus the actual dissipation is 
nonlocal and quite complicated. 

Now we consider the dissipative force on the ions. An exact analog of 
the calculation leading up to Eq. (42) gives the result that the 
dissipative force density on the ion fluid is given by 
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For the case of positive phase velocity, 1^/^ < 0 for a thermal ion 

distribution, so that the force is still a retarding force. Assuming 

w/k is roughly constant as a function of k, the perturbed ion pressure in 

the x domain has the same form as given in Eq. (45). Depending on the ion 

temperature, the corresponding to the ion Landau damping can be large or 

small. If the electron temperature is comparable to the ion temperature, 

the damping length is comparable to iL, implying that v < (—) ui 

u x>i v tn p6 

Hence for either electron or ion dissipation, values for are quite 
large, easily large enough that the electrostatic potential is highly 
localized in the overall magnetic structure. Furthermore, the dissipation 
is not anomalous, in the sense that it does not arise from any instability, 
as does the dissipation which generates the transverse v. 

Since ion Landau damping heats the ions, and electron Landau damping 
heats the electrons, a reasonable guess for is 
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where we have taken aj/k ■ V, the flow speed. 

We now continue this section by discussing other aspects of the Landau 
damping model. As complicated as say Eq. (44) for is, it is still a 
fairly coarse approximation to the actual nonlocal dissipation for several 
reasons. First of all, Eq. (44) is derived from linear theory, whereas the l 

shock wave is inherently nonlinear. Secondly, since it is nonlocal, 

perturbed velocities everywhere contribute to the perturbed dissipative 
pressure at each point. However, only a portion of the plasma is 
electrostatic; a large part of it is dominated by the magnetic structure, 
and this was not treated correctly at all in the calculation of P .. 

Third, the Fourier transform technique used to calculate P g works only 
if ? e is either periodic or square integrable. However, a shock profile is 
neither so it is not clear just how valid the Fourier representation is. 

For instance, if V g (x) in Eq. (44) has a typical shock structure, then the 
integral in Eq. (44) diverges for large x. Thus the linear theory in the 
Fourier domain gives an infinite perturbed pressure at large x, that is, 

the region of x furthest from the shock. It is the divergence which is at 

26 

the root of the claim by Pfirsch and Sudan that Landau damping cannot 
form a steady state shock. We feel that this divergence is clearly 
unphysical, so either nonlinear contributions play an important role in 
determining P g or else the Fourier transform scheme must be modified. 

Hence, since the actual dissipation mechanism is so complicated, a 
simple fluid approximation to it could be extremely useful. The only four 
possible conventional dissipation mechanisms are thermal conductivity, 
viscosity, transverse resistivity, and longitudinal resistivity. Thermal 











conduction can almost certainly be eliminated as a mechanism to form the 


shock because there does not seem to be any way that diffusing temperature 
can slow the forward motion and turn the upstream kinetic energy into 
downstream temperature. Viscosity can form the shock. However, there is 
no maximum Mach number for a viscous shock, no way that it is related to 
the breakdown of quasi-neutrality, and no way that it is related to ion 
reflection. Since all of these are known to be important in experiments, 
viscosity does not appear to be the best fluid description of the 
dissipation either. Transverse resistivity is not only a good 
approximation, in all likelihood it is the dissipation mechanism, since 
transverse electron-ion streaming can give rise to instabilities which 
exchange momentum between electrons and ions.^ However, as we have seen, 
transverse resistivity can only form the shock for Class B plasmas. 

This leaves longitudinal resistivity. While the actual dissipation 
mechanism is considerably more complex than this, it is a very good 
approximation to the physics. It is directly related to the breakdown of 
quasi-neutrality, it has an upper Mach number limit corresponding to ion 
reflection, it predicts a very narrow isomagnetic potential jump at the 
shock crest, and it has a lower Mach number limit corresponding to the 
maximum Mach number for transverse resistivity. All of these features are 
confirmed by the laboratory experiments we have cited. Furthermore, while 
longitudinal resistivity is simple, the two parameters which characterize 
it, and can be estimated from a coarse approximation to the actual 
dissipation mechanism. 

Let us conclude this section by discussing an important additional 
constraint on the physics. As we have seen, our basic model is that the 
flow steepens up until an ion acoustic wave is generated. The Landau 
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damping of this ion acoustic wave then provides the dissipation which forms 
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the shock. However, for the theory of Landau damping to be valid, the 
trapping time (for say the ions) 
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where $ is the electrostatic potential, must be longer than the correlation 
time for a particle with velocity v. 
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However, if the shock is simply a steady state potential profile, t =■ 00 

c 

for particles at the shock velocity. Hence the entire concept of Landau 
damping breaks down and the ion dynamics is dominated by trapping in a 
stationary potential. Thus in order for the Landau damping model to be 
valid, either the shock must have a temporally fluctuating structure, or 
else the total time of the experiment must be less than a trapping time. 
It is very likely that the shocks do have a fluctuating structure if only 
because in the broader part, dominated by transverse resistivity, the 
resistivity is probably caused by Debye length scale streaming 
instabilities. The spectrum of unstable waves will almost certainly 
overlap the electrostatic subshock and give rise to overall fluctuations in 
time. 
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V. Results 


In this section we present results of the calculations outlined in 
Section III. There are two classes of calculations. First, we consider 
parameters like those examined experimentally in Ref. 22. Here we model 
the longitudinal resistive heating as going entirely into the ions, as 
would be the case for strong ion Landau damping (that is not much less 
than T e ). Secondly, we consider a nitrogen plasma with parameters like 
that given in Ref. 21. More important, however, we consider no resistive 
ion heating. 

A. The Case of Strong Ion Heating 

In this case, we consider a shock in a hydrogen plasma with parameters 

as given in Ref. 22. That is, the upstream density is 4.6x10^ cm - ^ and 

the upstream field is 730G. The ion y is taken as equal to 5/3, and T g » 

T^ » 1 ev. As Is conventional in resistive shocks we assume a length of 

roughly 10 c/co so that cv/co V. * 10 where we use upstream values for the 
pe pe A 

plasma frequency and Alfven speed. Although it is often assumed that for 
transverse resistive shocks there is no resistive ion heating, it is not 
easy to envision a transverse streaming instability which gives none. 
Also, a small amount of resistive ion heating would be extremely difficult 
to measure, and it might in fact be very important in determining the 

Q 

resistive shock width. In addition, the ion temperature greatly affects 
the ion Landau damping rate in the electrostatic portion of the shock. 
Accordingly we take a * 0.2, so there is 5 times as much resistive electron 
heating as resistive ion heating in the transverse part of the shock. 










Since the transverse resistive shock always has a monatonic structure, 
the critical Mach number for transition to the Class C downstream state 
can be determined from the Rankine-Hugoniot conditions. For the 
parameters cited, this critical Mach number turns out to be 2.64. (Here 
the Mach number is defined as the shock velocity divided by the upstream 
magnetosonic speed). For Mach numbers below this, the shock is quasi¬ 
neutral and its structure is determined as specified in Section III.A. 

Above this Mach number, the quasi-neutral assumption breaks down. For 

larger Mach numbers, the shock structure is calculated as described in 

Secs. III.B and C. We assume, as seems consistent with experimental data, 

that in the electrostatic portion, ion Landau damping dominates. Thus we 

take o. * <*» and v (—. The upper critical Mach number, the 

£ £ 2 '-m' pe 

transition from Class C to Class D cannot be determined from Rankine- 
Hugoniot conditions alone for two reasons. First of all the electrostatic 
part of the shock generally is not monatonic, but overshoots the final 
density, and then decays down to the downstream value. Secondly, even if 
the shock were monatonic, the Rankine-Hugoniot condition gives only T e + 
T^, whereas the transition from Class C to Class D depends only on T^. 
Therefore this upper critical Mach number depends on how the resistive 
heating is partitioned between electrons and ions. For the parameters we 
have chosen, we find that this upper critical Mach number turns out to be 
about 4.46. This agrees reasonably well with the upper limit Mach number 
of 4.25 + 0.2 in Ref. 22. In Figs, la and lb are shown magnetic field and 
density profiles for a Class B shock with m ■ 2.5. Distances are given in 
centimeters. In Fig. 2a and 2b are shown magnetic field and density 
profiles in the transverse part of a supercritical shock having m * 4.4. 
The transverse part ends when the magnetic field achieves its downstream 



















value. At this point the electrostatic subshock is excited. The electron 
and ion density as a function of x in this region is shown in Fig. 2c. 
Again, distances are in centimeters. In all calculations, downstream 
values are carefully compared with Rankine-Hugoniot conditions as an 
independent check on the calculations. 

The actual ion heating downstream was measured in Ref. 22 by measuring 
n e (T g + T^) magnetically, and then measuring n g and T 0 by laser 

scattering. The results, along with our calculated results, are shown in 
Fig. 3. The two basically agree, although the experiment shows somewhat 
more ion heating than what is calculated. 

Another important experimental result from Ref. 23 is that as the Mach 

2 2 

number increases, the value of 2e$/M(V o - V^) decreases. Here (j> is the 
total potential jump and V Q and are upstream and downstream flow 

speeds. This goes hand in hand with the idea that as the Mach number 
increases, ions are slowed down not only by the electric field, but also by 
some additional dissipative effect. Our dissipative model, longitudinal 
resistivity also provides additional ion slowing down. In Fig. 4 are shown 
variations of this parameter with Mach number as measured by Eselevich et 
al., and as we calculated assuming upstream parameters in Ref. 22. Again 
the agreement is quite good. 

B. Shocks with No Ion Heating 

In this section we consider shocks with no ion heating, a * = 0. 

This corresponds to zero temperature ions, or else a waterbag ion 
distribution in one dimension. We consider the latter, so y^ “ 3. Since 
the electrons and ions have different values of y, the Rankine-Hugoniot 
conditions cannot be calculated in the conventional way. However, by 
making use of the fact that a * 0, and the ions only heat adiabatically, 
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the conservation relations can be calculated. This is done in Appendix 


C. For upstream fluid parameters, we use those in Ref. 21, a singly 
ionized nitrogen plasma, n - lO^" 5 , B = 2KG, T g = 2 ev, = 0.5 ev, 
cv/oip e V A ■ 10, and v^ =■ u)^, corresponding to electron-Landau damping. 
Then, Rankine-Hugoniot conditions show that the Mach number for transition 
from Class B to Class C is M = 2.68. There is no second critical Mach 
number for transition from Class C to Class D predited by Rankine-Hugoniot 
simply because there is not enough ion heating. However, the shock is not 
monotonic, but oscillatory, so that in an overshoot it is possible to reach 
the second critical point V = 3T^/M. We find this singular point is 
reached at m * 3.73. In Fig. 5a and 5b are shown magnetic field and 


density profiles versus position in the quasi^neutral part for a shock 
with m ■ 3.5. In Fig. 5c is the ion density versus x in the electrostatic 

part. 
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VI. 


Conclusions 


Our first conclusion is that the experimental evidence strongly 
indicates the existence of a significant intermediate range of Mach numbers 
between shocks formed by transverse resistivity and shocks formed by 
reflected ions. In this range, according to Ref. 23, the breakdown of 
quasi-neutrality plays an important role and an electrostatic subshock 
forms on the crest of the magnetic structure. Second, we have shown that 
the singular behavior of the quasi-neutral fluid equations is not related 
to ion reflection but rather to the breakdown of quasi-neutrality. In the 
intermediate regime then, an ion acoustic wave train is set up. We 
speculate that the Landau damping of this acoustic wave is the dissipation 
mechanism which forms the shock. Third, we have investigated longitudinal 
resistivity as a simple fluid model for the dissipation mechanism for shock 
formation. We conclude that it is a reasonable model. It is important 
only between the upper Mach number for transverse resistivity, and the 
lower Mach number for a fluid picture of ion reflection. It also is 
directly related to the breakdown of the quasi-neutrality assumption, and 
it predicts an electrostatic subshock at the crest of the magnetic 
structure. 
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Appendix A 


In this appendix, we show that the scheme for determining the shock 
structure, as specified in Section III.C satisfy the overall conservation 
relations. Three of the conservation relations, that for mass, momentum 
and energy are satisfied at each point in the shock. The fourth relation 
uB » U Q B 0 is only satisfied downstream. 

Consider now the conservation equations for momentum and energy 
flux. Each flux is the sum of a fluid and electromagnetic component. In 
the quasi-neutral region, the electrostatic part of the momentum flux is 
neglected, but the magnetic part of the momentum flux, B^/8 it and the 
electromagnetic part of the energy flux BE q are each included in the 
formulation. In the quasi-neutral region then, n, V, and T are determined, 
in part, so that particle, momentum and total energy flux are conserved. 
At the switch point, B takes on its downstream value and remains 
constant. Therefore In the electrostatic region, the magnetic part of the 
momentum flux and electromagnetic part of the energy flux are constant and 
have the proper downstream value. 

In the electrostatic region, the formulation leads to a conservation 

law for total energy flux and a conservation law for total fluid momentum 

E 2 

flux plus electrostatic momentum flux - . Downstream, of course, E =* 0 
and the fluid momentum and energy fluxes are as given by the Rankine- 
Hugoniot condition. 

Now consider the fluid energy flux as one moves back upstream toward 
the switch point in the electrostatic region. Since it is constant, it 
matches smoothly onto that in the quasi-neutral regime. Next consider the 
momentum flux. Since the formulation for the electrostatic region contains 







a momentum flux - , and the quasi-neutral formulation does not, the 


momentum fluxes do not quite match at the switch point if they have the 

E2 

same values downstream. However, if -5— is a very small part of the 

OTT 

momentum flux at the switch point, the error in momentum flux is very 

small. Using the expression for E at the switch point, Eq. (2), and 

considering say the ratio of electrostatic to thermal momentum flux, we 
E^ 2 2 2 

find -x- /nT ~ feu v/vc k_l << 1. In all of our calculations of shock 
8 it e 1 pe D' 

structure, we have checked that the electrostatic component of the momentum 
flux is negligible at the switch point, so that total momentum is conserved 
throughout the shock. Of course, while the electrostatic momentum flux is 
very small at the beginning and end of the electrostatic shock, it can be 


















In this appendix, we calculate Rankine-Hugoniot conditions for the < 

case of Y e * y , * 3, but where the ions are only heated adiabatically 

so that T i - T io (IN 2 - Then denoting upsteam quantities with a subscript 
1 o 

zero and downstream quantities with no subscript, the conservation I 

equations for momentum and energy flux reduce to two equations for n and 
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In writing Eqs. (Bl) and (B2) we have used the fact that n Q V 0 = nV and V Q B 0 
■ VB. Now each equation is linear in T g , so solving for T £ from one and 
inserting in the other, we find 
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where u * — and V; ■ B /4irn M. Factoring out the obvious root u * 1, we 
n Ao o o 
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arrive at a cubic equation for u 
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This equation also has a root u ■ 1 if MV * MV, + T T + 3T, as can be 

o Ao J eo 1 

easily verified. The procedure then is to solve Eq. (B4) numerically 

starting from this root and following it as V Q increases. Assuming a 

singly ionized nitrogen plasma with n Q - 10 15 , T g0 = 2ev, T lo = 0.5ev, B * 

2KG, Fig. B1 shows a plot of downstream values of density and electron 

2 5 

temperature versus magnetosonic Mach number m ■ V/fV. + -=- T /M + 
















Downstream electron temperature (dotted) and density ratio 
(solid) for a shock with = 5/3, * 3 and adiabatic ion 

heating. 
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